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Abstract 

Background: The shrimp Nematocarcinus lanceopes Bate, 1888 is found in the deep sea around Antarctica and sub- 
Antarctic islands. Previous studies on mitochondrial data and species distribution models provided evidence for a 
homogenous circum-Antarctic population of N. lanceopes. However, to analyze the fine-scale population genetic 
structure and to examine influences of abiotic environmental conditions on population composition and genetic 
diversity, a set of fast evolving nuclear microsatellite markers is required. 

Findings: We report the isolation and characterization of nine polymorphic microsatellite markers from the 
Antarctic deep-sea shrimp species Nematocarcinus lanceopes (Crustacea: Decapoda: Caridea). Microsatellite markers 
were screened in 55 individuals from different locations around the Antarctic continent. All markers were 
polymorphic with 9 to 25 alleles per locus. The observed heterozygosity ranged from 0.545 to 0.927 and the 
expected heterozygosity from 0.549 to 0.934. 

Conclusions: The reported markers provide a novel tool to study genetic structure and diversity in Nematocarcinus 
lanceopes populations in the Southern Ocean and monitor effects of ongoing climate change in the region on the 
populations inhabiting these. 
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Findings 

The shrimp Nematocarcinus lanceopes Bate, 1888 [1] is 
well-lcnown from the deep sea around Antarctica, the 
sub-Antarctic islands and the adjacent deep-sea basins 
of the Southern Ocean [2]. It has a wide vertical distri- 
bution from the Antarctic continental slope down to the 
Southern Ocean abyssal plains at depths of 4,000 m [3]. 
Previous work using mitochondrial and nuclear se- 
quence data provided clear evidence for a single homo- 
geneous circum-Antarctic N. lanceopes population [4]. 
These findings were also supported by species distribu- 
tion modeling data, showing a large connected habitat 
around the Antarctic continent and the sub-Antarctic 
islands [5]. Currently, the Antarctic continent faces in 
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some areas dramatic and most rapid effects of climate 
change, e.g. atmosphere and ocean warming at the 
Antarctic Peninsula and in part the sub-Antarctic 
regions [6]. However, other regions appear quite stable, 
e.g. eastern Antarctica and the surrounding deep-sea 
basins [6,7]. The increase in temperature leads to strong 
declines in winter sea ice extent in the Southern Atlantic 
and is regarded as a likely cause for the observed dra- 
matic population decline in Antarctic krill populations 
(see [8]). The decline of this ecological key species has 
severe impacts on other animals in the Southern Ocean 
ecosystem that rely on this food source [9]. The 
consequences of this drastically changing environment 
on most, in particular benthic deep-sea organisms, are 
largely unknown and difficult to assess by direct 
observations. However, the expected responses of species 
and communities to climate change are critically 
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important for conservation plans and therefore analyses 
of population processes such as gene flow, local adapta- 
tion and loss of genetic diversity are of particular interest 
(see the Convention on the Conservation of Antarctic 
Marine Living Resources). The reported microsatellite 
markers allow monitoring the genetic composition of 
populations, and identifying declines in genetic diversity. 
Further they enable the identification of genotypes with 
selective advantages in regions with increased warming 
or higher rates of primary production [10]. 

We extracted DNA from pleon muscle tissue using a 
commercial extraction kit (QIAmp® Blood and Tissue 
Kit, Qiagen GmbH, Hilden, Germany), following the 
manufacturer's protocol. Two different approaches of 
microsatellite isolation were performed: First we created 
a small microsatellite enriched library following the Re- 
porter Genome Protocol (RGP) [11]. As a result of this 
approach we have detected 39 putative microsatellite 
loci from 96 clone sequences. Out of the 39 tested pri- 
mer pairs 18 loci amplified successfully. After fragment 
analysis and testing for polymorphic sites, only one 
microsatellite marker could be retained due to a lack of 
polymorphism, stutter bands or unreliable amplification. 
Therefore, we then used data from a next generation 
454 database created in an earlier study to identify 
putative microsatellite loci in Nematocarcinus lanceopes 
(see [12] for information on library creation and 



microsatellite statistics). The library was sequenced on a 
GS-FLX sequencer (Roche, 1/8 of a 454 plate). Coverage 
estimates for raw reads ranged from 0.11-1.33% [12]. 
Sequence assembly was performed using the software 
Geneious 5.5.6. [13]. Microsatellites were detected with 
the PHOBOS program v3.3.10 integrated as a plugin in 
Geneious (www.ruhr-uni-bochum.de/ecoevo/cm/cm_phobos. 
htm). Suitable oligonucleotide primers were designed using 
Primer3 [14] and forward primers tailed with either T7 or Sp6 
primer sequence for economic fluorescent labeling of 
PGR fragments [15]. A "pigtail" sequence was added to 
problematic reverse primers to reduce stutter bands [16] 
(see Table 1). We tested 57 primer pairs in the 454 ap- 
proach whereof eight loci could be retained after testing 
for polymorphic sites. 

Amplification success and polymorphism of the 
reported microsatellite loci were assessed using 55 
individuals of N. lanceopes from sampling sites all 
around the Antarctic continent illustrated in Figure 1. 
PGR amplification of loci was performed with the 
Qiagen Multiplex PGR kit (Qiagen Inc., USA) in 10 ^1 
reaction volumes using a GeneAmp 2720 Thermo Cycler 
(Applied Biosystems, USA). The temperature profile for 
amplification was as follows: (1) initial denaturation at 
95''G for 15 min, (2) 20 cycles: 95''C for 60 s, 58''C for 
30 s and 72°C for 30 s followed by (3) 10 cycles: 95°C for 
60 s, at 53°C for 30 s and 72''C for 30 s and (4) final 



Table 1 Characteristics of nine polymorphic microsatellite loci for Nematocarcinus lanceopes tested on 55 individuals 



Locus 


Repeat motif 


Oligonucleotide primer sequence 5' 3' 


Fragment length (bp) 


Na 


Ho 


He 


^Null Alleles 


A06 


(AG)25 


Fw: GTCCTGAGTAATCGGCTCAGCTCT 
Rev: ACCCAGTTGGAAGCTGTTCTGAG 


1 87-233 


17 


0.927 


0.924 


0 


NL37 


(AAT)19 


Fw: GGGTrTAGGAGGAGrrrCGGGAC 
Rev: ACGAGrrGGAATGGGGCTGATG 


1 86-228 


12 


0782 


0.832 


0.017 


NL16 


(AG) 12 


Fw: TCAATTGTCCGGGACGCAAATGT 
Rev: AGTACATGGGCCACTAACTCCG 


150-182 


13 


0.545 


0.549 


0 


NL44 


(ACAG)II 


Fw: AATGGAGTGCAATGACGCTTGG 
Rev: TCGCAGI 1 IGI 1 1 1 AGAGGGAGC 


254-328 


12 


0.818 


0.881 


0.026 


NL51 


(AT)12 


Fw: TGATGACAGGGAHTGTCTTTCG 
Rev: TCCCCATTTGTACGCTATCC 


144-165 


9 


0.727 


0.824 


0.052 


NL53 


(AC) 11 


Fw: ACAGTACACAGGCTACATAC 
Rev: ATCTTCATGTTATGCCCTCrAG 


149-179 


14 


0.800 


0.796 


0.001 


NL55 


(AG) 11 


Fw: ACGCGAACAGTGCTAAGAAGAC 
Rev: CACCACAGCAAGGAACCTCC 


162-216 


25 


0.927 


0.940 


0 


NL56 


(AG) 12 


Fw: AGTGAAAAGACTCAAATTCCTTGG 
Rev: CATTACTGCTTCCTTCCTCTC 


151-201 


23 


0.764* 


0.934 


0.077 


NL49 


(AG) 11 


Fw: ACTCTACirrGGCTTTCTCCCTC 
Rev: ACACGGGTCTTCCTGAGTGTTG 


155-183 


11 


0.927 


0.861 


0 



Nfi number of alleles; Hq observed heterozygosity; expected heterozygosity; F,vu// Alleles frequency of null alleles. 

* Indicates a locus that deviates from HWE after Bonferroni correction (9 comparisons; significance threshold = 0.005}. Forward primers of the loci A06, NL44 and 
NL55 were labelled with a T7 tail (5'-TAATACGACTCACTATAG-3') all other forward primers were labelled with a Sp6 tail (5'-GATTTAGGTGACACTAT-3'). A "pigtail" 
(5'-GTTTCTT-3') was added to the reverse primers of the loci NL49, NL51, NL53, NL55 and NL56 to reduce stutter bands [16]. 
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extension at 72°C for 40 min. We used slightly modified 
protocols for steps 2 and 3 for amplifying three micro- 
satellite loci: the loci NL55 and NL56 were amplified 
with annealing temperatures of 60°C (20 cycles) and 
55°C (10 cycles) and locus NL44 with annealing 
temperatures of 55°C (20 cycles) and 53°C (10 cycles). 
PGR products were multiplexed and added to a mixture 
of Hi-Di™ formamide and size standard (Gene Scan 
500 LIZ) and run on a 3500XL Genetic Analyzer capil- 
lary sequencer (Applied Biosystems, USA). To identify 
the obtained alleles based on the resulting electropherograms, 
the software GENEMAPPER v4.0 (Applied Biosystems, 
USA) was used. The software INEst [17] was used to esti- 
mate allele frequencies, observed and expected heterozy- 
gosity of each locus, and to assess the proportion of 
possible non-amplif}dng ("null") alleles. Of the 96 tested 
primer pairs, nine loci reliably amplified the desired locus 
and were polymorphic with 9-25 alleles per locus. The 
observed heterozygosity (Hq) ranged from 0.545 to 0.927 
and the expected heterozygosity (He) from 0.529 to 0.934. 
Null allele frequencies ranged between 0 and 0.077. The 
software GENEPOP v4.1.4 [18] was used to test for 
deviations from Hardy- Weinberg equilibrium (HWE) and 
for linkage disequilibrium. After Bonferroni correction for 
multiple comparisons, one locus showed significant 
deviation from expectations under HWE (see Table 1). 
Linkage disequilibrium was not detected for any of 
the pairwise loci comparisons. Nucleotide sequence 



information of the described loci is deposited in the Dryad 
Repository (http://datadryad.org/) and accessible via the 
link http://dx.doi.org/10.5061/dryad.gg3rn. 

The reported set of polymorphic microsatellite loci from 
Nematocarcinus lanceopes will help to analyze the fine-scale 
genetic structure around the Antarctic continent. These 
results will help us to test, whether there is regional differen- 
tiation that we have not yet detected relying on a single mito- 
chondrial marker only [4], whether gene flow is maintained 
primarily by the Antarctic Circumpolar Current, whether 
populations in regions under severe climate change (e.g., the 
West Antarctica) experience population declines or possibly 
growth, and whether patterns of local adaptations can be 
detected in the different regions. The markers will also allow 
monitoring the effects of the rapid climate change on the 
population genetic structure of this species. 
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